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ABSTRACT 

Context. We present the analysis of the galaxy structural parameters from Ha3, an Ho- narrow-band imaging follow-up survey of 
~ 800 galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) in the Local Supercluster, including the Virgo 
cluster, and in the Coma Supercluster. 

Aims. Taking advantage of Ho-3, which provides the complete census of the recent star-forming, Hi-rich galaxies in the local universe, 
we aim to investigate the structural parameters of the young (< 10 Myr) and the old (> 1 Gyr) stellar populations. By comparing the 
sizes of these stellar components, we investigated the spatial scale on which galaxies are growing at the present cosmological epoch 
and the role of the environment in quenching the star-formation activity. 

Methods. We computed the concentration, asymmetry, and dumpiness (CAS) structural parameters for recently born and old stars. To 
quantify the sizes we computed half-light radii and a new parameter dubbed EW/r based on the half-light radius of the Ha equivalent 
width map. To highlight the environmental perturbation, we adopt an updated calibration of the Hi-deficiency parameter (Def H] ) that 
we use to divide the sample in unperturbed galaxies (Def H i < 0.3) and perturbed galaxies (Def HI > 0.3). 

Results. The concentration index computed in the r band depends on the stellar mass and on the Hubble type these variables are 
related because most massive galaxies are bulge dominated therefore highly concentrated. Going toward later spirals and irregulars 
the concentration index and the mass decrease along with the bulge-to-disk ratio. Blue compact dwarfs (BCDs) are an exception 
because they have similar mass, but they are more concentrated than dwarf irregulars. The asymmetry and the dumpiness increase 
along the spiral sequence up to Sc-Sd, but they decrease going in the dwarf regime, where the light distribution is smooth and more 
symmetric. When measured on Ha images, the CAS parameters show no obvious correlations with Hubble type. Irrespective of 
whether we used the ratio between effective radii or the EW/r parameter, we found that the concentration index is the main parameter 
that describes the current growth of isolated galaxies but, for a fixed concentration, the stellar mass plays a second-order role. At the 
present epoch, massive galaxies are growing inside-out, conversely, the dwarfs are growing on the scale of their already assembled 
mass. 

Key words. Galaxies: clusters: individual: Virgo - Galaxies: clusters: individual: Coma - Galaxies: fundamental parameters structure, 
sizes - Galaxies: star formation 



* Observations taken at the observatory of San Pedro Martir 
(Baja California, Mexico), belonging to the Mexican Observatorio 
Astronomico Nacional. 

** Tables Al and A2 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+ A/ 
*** Hubble Fellow 



1. Introduction 

Understanding the physical structure of galaxies is a key to shed 
light on the assembly of the baryonic matter in the Universe 
during the cosmic growth of large-scale structures. The first 
and most widely used way to describe the physical structure 
of galaxies is the morphological classification. Since the 1920s 
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several morphological schemes have been proposed and adopted 
(Hubble 1936; de Vaucouleurs 1959). Although the adopted cri- 
teria were somewhat subjective, noticeably the different classi- 
fications proposed by the most distinguished astronomers were 
consistent. In recent years, the enormous size of galaxy samples 
(e.g., those produced by the Sloan Digital Sky Survey, SDSS) 
made the visual estimate of the morphology virtually impossi- 
ble (with the exception of the Galaxy Zoo project by Lintott et 
al. (2008), which aims at a human classification using a quar- 
ter of million of volunteers). Nevertheless, the improvements of 
computer science made possible the advent of automated quanti- 
tative morphological classification methods, which are trying to 
make the scheme proposed by Hubble "quantitative" using ob- 
jective parameters (e.g. Cheng et al. 201 1). One alternative is to 
describe a galaxy parametrically, by modeling the distribution of 
light as projected onto the plane of the sky with a prescribed suite 
of analytic functions. Following this parametric modeling, the 
light profile can be decomposed into a bulge and a disk that can 
be fitted with a Sersic function (Sersic 1968) and an exponen- 
tial, respectively. However, Laurikainen et al. (2005) and Peng 
et al. (2010) have shown that simple bulge + disk models are 
inadequate for dwarf and irregular galaxies. 

More recently, a new non-parametric approach to galaxy 
morphological classification is emerging. This scheme does not 
assume a particular analytic function for the galaxy's light distri- 
bution and therefore may apply to irregulars as well as to galax- 
ies on the Hubble sequence (Lotz et al. 2004). In addition, these 
methods are less computationally intensive than profile fitting 
and can be used even at high redshift where the resolution and 
signal-to-noise are too low for a parametric (or even visual) clas- 
sification. One of the most popular schemes is the concentra- 
tion, asymmetry, and dumpiness (CAS) system developed by 
Conselice (2003), inspired by other seminal works (Abraham 
et al. 1996; Takamiya 1999; Bershady et al. 2000; Conselice et 
al 2000). Although several works use the CAS classification at 
high redshift (e.g. Conselice et al. 2003, 2007, 2008; Cassata et 
al. 2005), systematic attempts to quantify the CAS parameters 
at low redshift are rare in the literature. Previous efforts have 
been carried out by Hernandez-Toledo et al. (2008), applied to 
isolated galaxies selected from the CIG catalog (Karachentseva 
1973). The first goal of the present work is to produce a reference 
for the CAS parameters computed on a well-selected, complete 
sample that is representative of late-type galaxies at z ~ 0. 

Since Walter Baade (1944) discovered that young Pop-I stars 
have a broader spatial distribution than older Pop-II in spiral 
nebulae (made possible by wartime blackouts, which reduced 
the light pollution at Mount Wilson), it became increasingly 
clear that galaxies grow inside-out (see e.g. Munoz-Mateos et al. 
2007). Owing to the extensive sky coverage of the SDSS survey 
(Abazajian et al. 2009) and of the GALEX (Martin et al. 2005) 
mission, even an eyeball inspection of pseudo-color images of 
galaxies, especially the most inclined ones in the plane of the 
sky, reveal bluer colors in the outer wings, suggesting younger 
stellar populations (see e.g. de Jong 1996). However, a caveat 
is that radial gradients of metallicity (Skillman et al. 1996) and 
dust extinction would eventually artificially enhance any color 
gradient induced by a pure change in stellar age. Resolved stellar 
populations have been studied with the Hubble space telescope 
(HST) and revealed, for example, the inside-out growth of the 
disk galaxy M33 (Williams et al. 2009). A definitive proof of 
inside-out build-up of galaxies would need to be derived from 
a direct comparison of the radial extent of old stars (traced for 
example by r-band observations or redder) with that of stars as 
young as 10 6 ~ 7 yr, such as those revealed by Ha observations. 



Ground-based integral-field spectroscopy has made this possi- 
ble and demonstrated the power of spatially resolved emission 
line surveys at low- and high-redshift (e.g. Forster Schreiber et 
al. 2009; Law et al. 2009; Sanchez et al. 2012), but it has been 
limited to small samples due to the required extensive obser- 
vational effort. Recently, a comparison between the r- and Ha 
scale-lengths of galaxies at z — 1 was carried out by Nelson et 
al. (2012), using spatially resolved HST grism spectra from the 
3D-HST survey (Brammeret al. 2012), showing that the effect is 
mild at best. In addition to the CAS parameters, we present here 
the size measurements of local star-forming galaxies as traced 
by different stellar populations. 

The outline of the paper is as follows: Section [2] introduces 
the selection criteria and the completeness of the samples used 
in this work. In Section [3] the CAS parameters are reviewed, 
together with a description of computation details and validation 
tests for the specific software we developed for this paper. In 
Section |4] the methods used for the size measurements of the 
spatial extent of galaxies are presented. Finally, the discussion 
of our results follows in Section Throughout the paper we 
adopt a flat ACDM cosmology, with Ho = 73 km s _1 Mpc~' and 
Q.A = 0.7. All magnitudes are given in the AB system unless 
explicitly noted otherwise. 

2. Samples 

2.1. Local Supercluster 

The first sample analyzed in this work is drawn from the 900- 
square-degrees region between ll h < R.A. < 16 h and 4° < 
Dec. < 16°, which covers the Local Supercluster, including the 
Virgo cluster. This region has been fully mapped by ALFALFA 
(Giovanelli et al. 2005), which provides us with a complete sam- 
ple of Hl-selected galaxies with masses as low as 10 77 M Q 
(Haynes et al. 2011). 

Ha3 is a follow-up survey consisting of Ha-imaging obser- 
vations of ALFALFA detections with high signal-to-noise (typ- 
ically S IN > 6.5) and a good match of two independent polar- 
izations (code = 1 sources; Giovanelli et al. 2005, Haynes et al. 
2011). Ha observations of members of the Virgo cluster were 
presented in Gavazzi et al. (2002a,b, 2006), Boselli & Gavazzi 
(2002), and Boselli et al. (2002a,b). Images and fluxes are also 
publicly available via the GOLDMine web server (Gavazzi et al. 
2003). Ha data for galaxies in low-density regions of the Local 
Supercluster were presented in Gavazzi et al. (2012 hereafter 
Paper I). Outside the Virgo cluster the redshift window^ of Ha3 
is 350 < cz < 2000 km s , and we limit our selection to objects 
with HI fluxes Fm > 0.7 Jy km s , while in the Virgo cluster 
the velocity interval is extended to 350 < cz < 3000 km s , 
without any constraint on the HI flux, to map the cluster in its 
full extent. At the distance of 17 Mpc assumed for the Virgo 
subcluster A (Gavazzi et al. 1999), the adopted flux limit corre- 
sponds to an HI mass Mui = 10 7 7 M , while Ha3 is complete 
to M H i > 10 8 M Q (see Paper I). ALFALFA detects about 50% 
of the late-type galaxies in the Virgo cluster, either because the 
remaining are dwarf galaxies with an HI content below the lim- 
iting sensitivity or because they are highly Hi-deficient galaxies. 
Fig.Q~]illustrates the sky region covered by Ha3, which contains 
409 galaxies. Among these, 372 objects are late-type galaxies 
and are shown as red dots. 



1 The lower velocity limit is set to avoid galaxies whose Ha line falls 
on the steep shoulder of the transmission curve of the filter used at SPM 
(see Paper I). 
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R.A. (J2000) 

Fig. 1. Sky distribution (R.A. and Dec. not in scale) of the Local Supercluster galaxies with superposed the X-ray contours 
(Bohringer et al. 1994) highlighting the position of the Virgo cluster. Red points refer to the 372 late-type galaxies observed by 
Ha3 in the ALFALFA footprint. The 13 optically selected VCC galaxies followed-up in Ha are shown as blue points. 



In addition, 13 galaxies from the Virgo Cluster Catalog 
(VCC, Binggeli et al. 1985) that are undetected by ALFALFA 
have been previously observed in Ha. They are well-known 
highly Hi-deficient objects from previous deep pointed obser- 
vations taken at Arecibo by Helou et al. (1984), Haynes & 
Giovanelli (1986), Hoffman et al. (1987, 1989), and Gavazzi et 
al. (2005). Although these galaxies are not selected according to 
any criterion of completeness, we included them to fully exploit 
the dependence of the structural parameters on the amount of 
HI. These galaxies are shown in Fig. [T]as blue dots. The total 
number of late-type galaxies in this sample is 385. 

2.2. Coma Supercluster 

This second sample includes ~ 650 nearby galaxies, selected 
from ALFALFA in the "spring" sky of the Coma Supercluster. 
The traditional area of the Coma Supercluster (11.5 h < R.A. < 
13.5 h ; 18° < Dec. < 32°; 3900 < cz < 9500, hereafter re- 
gion 1) is a very complex environment composed of two ma- 
jor clusters, Coma and A 1367, and smaller groups and filaments 
that connect the clusters and build the Great Wall (Gavazzi et 
al. 2010). Only a small strip (hereafter region 2) between 10 h < 
R.A. < 16\ 24° < Dec. < 28°, and 3900 < cz < 950010ns- 1 
is currently included in the ALFALFA catalog (Haynes et al. 
2011). At present, half of the Coma cluster lies in the footprint 
of ALFALFA. We found 683 galaxies with high S/N and a good 
match of two independent polarizations. These are the same cri- 
teria as for the Local Supercluster selection but, because of the 
distance of Coma (lOOMpc) we applied no additional threshold 
on the limiting HI flux. The limiting sensitivity of ALFALFA 
at the distance of Coma is Mui ii m = 10 91 M o . Furthermore, the 
Ha3 survey is still not complete in this strip, where only 261 
galaxies have been followed-up in Ha at present. 

Because neither ALFALFA nor Ha3 are currently complete 
in the whole region under study, we complemented the HI se- 
lection with galaxies optically selected from the CGCG catalog 
(Zwicky et al. 1961-68) that were previously observed in Ha. 
We found 372 bright (m\i m = 15.7 mag veg!i ) late-type galaxies in 
the CGCG in region 1. Of these, 82 are in the region covered 
by ALFALFA and are Hi-detected therefore we refer to them as 



radio-selected galaxies. The remaining 290 galaxies lie beyond 
the ALFALFA footprint or are Hi-undetected by ALFALFA, 
therefore they are considered to be optically selected. The ma- 
jority of these optically selected galaxies (230) are detected in 
HI from previously deep pointed observations (see references in 
Sec 12. U . We recall that because the HI observations have been 
carried out with different integration times, this sample is not 
complete to a well defined Hi-mass limit. Nonetheless, the opti- 
cally selected galaxies have HI properties that mimick those of 
the radio-detected ones at the distance of Coma, and there is ev- 
ery expectation that most of them (~ 80%) will be detected by 
ALFALFA. The remaining (~ 20%) are Hi-deficient galaxies, 
for which ALFALFA is too shallow. Our combined sample is 
at least as deep as ALFALFA for the HI content of the galaxies. 
We refer to the detailed comparison of the optical vs. radio selec- 
tion performed by Gavazzi et al. (2013b, hereafter Paper III) for 
the other properties. According to these authors, our selection of 
bright objects from the CGCG catalog allowed us to build a sam- 
ple whose stellar masses, colors and star formation rates are not 
significantly different from the late-type galaxies selected from 
ALFALFA. In conclusion, we expect to introduce no biases by 
combining the two selection criteria. 

The Ha images in the Coma sample are taken from Ha3 
(Gavazzi et al. in prep., hereafter Paper V) (see Fig. |2]red points) 
or from previous observations of 155 optically selected CGCG 
galaxies by Gavazzi et al. (1999, 2002a, 2002b, 2006), Boselli & 
Gavazzi (2002), Iglesias-Paramo et al. (2002), and Cortese et al. 
(2006) (see Fig.|2]blue points). For the CGCG, images and fluxes 
are publicly available via the GOLDMine web server (Gavazzi 
et al. 2003). The number of galaxies (with available Ha images) 
in the Coma Supercluster sample is 413, which brings the total 
number of galaxies analyzed in this work to 798. 

Following the prescriptions adopted in the previous papers 
of the series, the stellar mass was derived from the SDSS i mag- 
nitudes and g - i color using the transformation 

MS = _L94 + °' 59 ' iS ~ + L15 ' Mi;)' (1) 

where logL,- is the /-band luminosity in solar units (logL,- = 
(Imag - 4.56)/ - 2.5). This is a modification of the Bell et al. 
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R.A. (J2000) 



Fig. 2. Sky distribution (R.A. and Dec. not in scale) of the Coma Supercluster galaxies. Black dots represent the optically selected 
galaxies (all types) selected from the SDSS database (ll' ! 30 m < R.A. < 13 /! 30 m , region 1) or from ALFALFA in the strip between 
24° < Dec. < 28° beyond the previous R.A. window (region 2). Red dots refer to the 261 galaxies observed in Ha3 in the ALFALFA 
strip. The 139 optically selected CGCG galaxies followed-up in Ha outside the ALFALFA strip and 16 in the ALFALFA strip are 
shown as blue dots. 



(2003) formula introduced for consistency with the mass deter- 
mination of MPA-JH ifl Morphological types are taken from the 
VCC and CGCG catalogs when available. For all the others the 
type was assigned after visual inspection of the SDSS pseudo- 
color images. Following a procedure similar to the one described 
in Gavazzi et al. (2010), giant late-type galaxies (Sa - Sd) were 
classified in steps of half Hubble class on a purely morpholog- 
ical basis. Among dwarf galaxies, blue compact objects (BCD) 
were classified according to their amorphous compact morphol- 
ogy and blue color, on the other hand, dlrr are blue but low sur- 
face brightness objects. 

3. The CAS System 

One of the most popular non-parametric sets of morphologi- 
cal measurements is the CAS system introduced by Conselice 
(2003). It is based on three indices: concentration C, asymmetry 
A, and dumpiness S . 

3.1. Computing CAS parameters 

Before structural parameters measurements were obtained, ba- 
sic digital image-processing was applied to all frames (i.e., 
bias-subtraction, flat-fielding, cosmic-ray removal, and sky- 
subtraction; see Paper I). In addition, any objects in the galaxy's 
image that are not part of the galaxy, such as foreground stars 
or foreground and background galaxies, were removed using the 
IRAF task imedit. CAS parameters in the Gunn r-band images 
were computed as described in the following subsections for all 
galaxies except for close pairs or those that are affected by haloes 
of bright stars. Furthermore, we disregarded galaxies whose sur- 
face brightness is so low that the centering and noise-correction 
routines never converge. These two problems prevented us from 
measuring structural parameters for 56 galaxies in our samples. 

The CAS parameters were also computed in the Ha images. 
The center of the images was set to be the same computed in the 

2 http://www.mpa-garching.mpg.de/SDSS/DR7/ 



r band since Ha is often centrally depressed (Shapley 201 1). In 
addition, the concentration index is affected if there is an active 
galactic nucleus (AGN), because in these cases the central emis- 
sion in the Ha images does not trace recent star formation but 
the ionization of the gas caused by the AGN. To overcome this 
problem, we used the nuclear spectra from Gavazzi et al. (201 1) 
to identify galaxies with an AGN. Then, the central region was 
masked in these images to exclude the effects of the AGN emis- 
sion. The standard radius for this mask was set to be 10 arcsec 
for nearby galaxies in the Local Supercluster and 2.5 arcsec for 
galaxies belonging to the Coma Superclustefl 

We also excluded from the analysis of the CAS parameters 
107 galaxies that are either undetected in Ha or that show only 
nuclear emission associated to an AGN. In the end, structural pa- 
rameters were measured on r-band images for 742/798 galaxies 
and on Ha images for 635/798 galaxies. 

3.2. Concentration index 

Concentration is defined in slightly different ways by different 
authors, but it basically measures the ratio of light within a cir- 
cular or elliptical inner aperture to the light within an outer aper- 
ture. The CAS system adopts the Bershady et al. (2000) defini- 
tion as the ratio of the circular radii containing 20% and 80% of 
the "total flux": 

C = 5x/ g(— I (2) 

where rgo and are the circular apertures containing 80% and 
20% of the total flux. As a definition of the total flux, Conselice 

3 The Coma Supercluster is five times farther away than Virgo, 
therefore the masking radius should be 2 arcsec. We used 2.5 arcsec 
(~ 1.5 x < PSF FWHM >) to take into account the average seeing con- 
ditions of our images. For extended low ionization nuclear emission re- 
gions (see Keel 1983; Hameed & Devereux 1999) which, in most cases 
are ionized by post-asymptotic giant branch (PAGB) stars (Stasihska et 
al. 2008), we used a masking radius twice as large, i.e., 20 arcsec at the 
distance of Virgo and 5 arcsec at the distance of Coma. 
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(2003) used the flux contained within 1 .5 times the Petrosian ra- 
diufl Tp (/so and r^o are considered only if they are larger than 
1.5 times the seeing, as suggested by Arribas et al. (2012), oth- 
erwise C is not given). For the measurement of C the galaxy's 
center coincides with the center that minimizes the asymmetry 
(see below). 

Yamauchi et al. (2005) pointed out that the concentration 
index can be overestimated for edge-on galaxies when circu- 
lar apertures are used. We performed multi-aperture photometry 
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Fig. 3. Difference between the concentration index in the r- 
band computed in elliptical and circular apertures as a function 
of £ = log 10 (a/Z?) (de Vaucouleurs et al. 1991) for the Local 
Supercluster (left) and the Coma Supercluster (right) samples. 



both in circular and elliptical apertures and we found a strong 
correlation between the axis ratio and the difference between the 
concentration index computed in circular (C) and elliptical (C e ) 
apertures, as showed in Fig. [3] In other words, edge-on spirals 
have an abnormally high concentration of light if defined with 
circular apertures, and the morphological classifications based 
on C will suffer from significant contamination. Using elliptical 
apertures remedies this problem, therefore our CAS results are 
based on the concentration index C e . 



3.3. Asymmetry index 

The asymmetry parameter A quantifies the degree to which the 
light of a galaxy is rotationally symmetric. A is defined as the 
difference in flux between the original image and the one rotated 
by 180° from its center, normalized to the flux in the original 
image (Abraham et al. 1996; Conselice et al. 2000): 



A 



I,ijm,J)-Rm(i,J)\ 



(3) 



where / is the original image and Ri&o is the image rotated by 
180°. A is summed over all pixels (/,/') within 1.5 r p of the 
galaxy center. 

The computation of the asymmetry is not robust with respect 
to the assumed photometric center, since a positional difference 
of 1 arcsec produces substantially different asymmetry values. 



The Petrosian radius is the radius r p at which the ratio of the surface 
brightness at r p to the mean surface brightness within r p is equal to a 
fixed value, typically set to 0.2 (Petrosian 1976). 



To overcome this centering problem, we defined the rotation cen- 
ter as the position that yields a minimum value of the asymme- 
try parameter. It should be noted that this issue especially affects 
the asymmetry parameter but, following Conselice (2003), Lotz 
et al. (2004), and Vikram et al. (2010), we used this centering 
method also compute the two other CAS parameters. 

Operationally, we made an initial guess starting from the 
photometric center and then iteratively computed the asymme- 
try in a 5 x 5 pixel grid. Unlike those of other authors (Conselice 
2003; Vikram et al. 2010), our grid steps are integer pixels Cotini 
et al. (2013) claim that the use of fractional pixels introduces a 
non-negligible smoothing that affects the CAS parameters. The 
Petrosian radius was recomputed after each step, until the asym- 
metry reached a minimum. The asymmetry computed in this way 
can be overestimated if the background contains structure (i.e., 
haloes of bright stars or poor flat fielding). These effects were 
corrected for by simultaneously performing the asymmetry mea- 
surement on both the source and a neighboring blank area of the 
image. Taking into account these noise corrections, the equation 
used to compute the asymmetry is 

, XijW,J)-Riw(i,j)\ Z U \N(iJ)-NR m (iJ)\ 

A = — - — -T— — -T— — , (4) 



Z;,;|/(U)I 



where N is the background image and NR is the rotated back- 
ground image. 

3.4. dumpiness index 

The dumpiness parameter S has been developed by Conselice 
(2003), inspired by the work of Takamiya (1999), to quantify 
the amount of light in small-scale structures. The dumpiness is 
defined as the ratio of the amount of light contained in high- 
frequency structures to the total amount of light in the galaxy. 
For ellipticals this ratio should be, and generally is, near zero. 

First, the original galaxy effective resolution was reduced 
to create a blurred image (B). This was done by smoothing the 
galaxy with a boxcar filter of width <x. Although in principle any 
scale can be used to measure various frequency components in a 
galaxy, we followed Conselice (2003) and Lotz et al. (2004) set- 
ting a smoothing filter size cr(I) = 1 /5 X 1 .5 X r p = 0.3 X r p . The 
effect of the blurring is to create an image whose high-frequency 
structures have been washed out. If a galaxy only consists of 
concentrated but structureless light (like a BCD, for example) 
it has a very low dumpiness. The blurred image was then sub- 
tracted from the original image (/). The dumpiness parameter S 
is defined as the flux in the residual image normalized to the flux 
in the original image. 

The central pixels within a circular aperture equal to the 
smoothing length 0.3 r p were excluded from the sum because 
the central parts of the galaxies are often unresolved and con- 
tain significant high-frequency structure due to finite sampling 
(Conselice 2003). Following the definition of Conselice, we also 
forced to zero any negative pixel in the residual map before com- 
puting S . This was achieved by substituting the original image 
with the one (I max ) obtained taking the highest value between / 
and B for each pixel. 

Even the dumpiness can be overestimated by noise effects 
caused by the graininess of the background. We corrected this 
source of noise in the same way as for the asymmetry, and the 
final equation used to compute dumpiness becomes 

c , n ZytWi, j) ~ B(i, J)] n Z Lj [N m . dx (i, j) - NB(i, j)] 
S = 10 — - - — 10 — - - — ,(5) 



Z/j[^max0'-;)] 
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where N mm is the maximum background image and NB is the 
blurred background image. 

3.5. Comparison with the literature 

The code developed for the morphological analysis was tested 
on 18 Local Supercluster galaxies in common with Frei et al. 
(1996) that have also been analyzed by Conselice (2003). These 
are all giant, nearby NGC galaxies with Mb < -20 mag veg - d . 

Our estimates of the CAS parameters for this subsample are 
listed in Table Q] along with details on the telescope, the pixel 
scale and the seeing conditions (Full width half maximum of the 
point spread function, PSF FWHM) of the images. The Frei et 
al. images have been taken with similar r filters, but in slightly 
poorer seeing conditions (PSF FWHM ~ 2" - 5"). 

The comparison between the three CAS parameters com- 
puted using our in-house codes and those given in Conselice 
(2003) are shown in Fig. [4] The discrepant red point in the A 
and S comparisons is VCC1987 (NGC4654), a giant Sc galaxy 
in the Virgo cluster. This galaxy appears to be far more asym- 
metric and clumpy in our analysis than in Conselice (2003). We 
have several reasons to believe that this galaxy is in fact asym- 
metric. First, in the analysis of Conselice, this galaxy seems to 
be less asymmetric and clumpy than the remaining Sc galaxies in 
the sample. Furthermore, it is known that this galaxy is strongly 
asymmetric in the stellar disk, in atomic gas (HI) (Chung et al. 
2007), in molecular gas (H 2 ) (Cayatte et al. 1990; Chung et al. 
2009), and also in the Ha velocity field (Chemin et al. 2006). 
For these reasons we excluded this galaxy from the comparison. 

The agreement for the other points is remarkable. We found 
that the average deviation between our values and the pub- 
lished ones for concentration, asymmetry, and dumpiness are 
-0.02 ± 0.18, 0.03 ± 0.03, and 0.06 ± 0.09. Similar dispersions 
were seen in previous comparisons. The dispersion between the 
CAS parameters reported by Vikram et al. (2010) from that of 
Conselice (2003) are -0.11+ 0. 14, 0.00 ± 0.04, and 0.06 + 0.09 
using the whole Frei et al. (1996) sample. The agreement in these 
results gives an idea of the robustness of our software, but also of 
the stability of the CAS parameters when using different source 
images. We verified that convolving our images to a common 
PSF (with FWHM = 3") does not significantly change the com- 
parison with published values. Any resolution effect, if present, 
is negligible due to the small seeing differences between our im- 
ages and those by Frei et al. (1996) and to the limited statistics. 
Only by comparing the results from our whole sample with those 
of Conselice (2003) a small systematic difference can be seen 
and ascribed to a resolution effect (see Sect. 15. II ). Indeed, Lotz et 
al. (2004) found that C, A, and S parameters are reliable within 
+ 0.3, 0.1, and 0.2 when the S/N per pixel, within 1.5 x r p , is 
<S/N> > 5. 



4. Size measurements 

To investigate on which spatial scale galaxies are growing, we 
parameterized the spatial distribution of the old and the young 
stellar populations. The stellar mass of a galaxy was traced us- 
ing the Gunn r band and the star-formation was traced by the 
Ha line emission. Some uncertainties that affect the size mea- 
surements are to be considered. First, Ha /16563 and the [Nil] 
doublet AA 6548, 6584 are not spectrally resolved by means of 
narrow band interferential filters. We do know the line ratios 
for a large number of objects from integrated drift-scan spec- 
troscopy (Gavazzi et al. 2004; Boselli et al 2012), although they 



are of limited use for our purpose, since line ratios almost cer- 
tainly vary with the position in the galaxy (Genzel et al. 2008; 
Yuan et al. 201 1). Therefore we did not apply any mean correc- 
tion. We tried to reduce this problem by masking the AGN emis- 
sion (dominated by [Nil]), where present, using nuclear spectra 
(see below). A second systematic uncertainty is dust attenuation. 
Differential extinction could affect the spatial distribution of the 
Ha emission (Boissier et al. 2004, 2007), and thus the size mea- 
surements. Finally, the Gunn r filter contains the Ha (and [Nil]) 
lines, which introduces a correlation between the two images. 
Nevertheless ,this effect is small because the lines' contribution 
to the stellar continuum images is lower than 10% of the total 
flux (Paper I). 

4.1. Effective radii 

The effective radius (or half-light radius) is a parameterization of 
a galaxy's light distribution. It is defined as the radius containing 
half of the total light in a specific band. The standard method for 
determining half-light radii of galaxies is to fit two-dimensional 
models to the galaxy light profiles (e.g. using GALFIT, Peng 
et al. 2002). However, this parametric method does not provide 
meaningful results for the Ha images because they are typi- 
cally clumpy and often have asymmetric or centrally depressed 
light profiles (Shapley 2011). For this reason we followed the 
method suggested by Forster Schreiber et al. (2009), Arribas et 
al. (2012), and Nelson et al. (2012). We measured effective radii 
for both the Ha images and the optical r-band continuum us- 
ing integrated curves of growth (CoG). Using this method, the 
effective radii depend on the limiting sensitivity of the images. 
However, we compared the effective radii measured from CoGs 
with those obtained fitting the radial light profiles by Cortese 
et al. (2012) on a sample not presented here (Fossati et al. in 
prep.), and found only small differences (~ 10%) between the 
two methods. Another uncertainty can arise if the limiting sen- 
sitivity of narrow-band Ha images and that of the r broad-band 
are different. The sensitivities in those two bands are compara- 
ble if the exposure times are in the same ratio as the filter band- 
passes. We have a range of these exposure time ratios spanning 
from 60% to 100%. We verified that using the galaxies observed 
in the correct ratio of the exposure times and simulating 40% 
shorter exposures in the narrow -band filter, the derived effective 
radius r e (Ha) typically varies by only ~ 5%. We conclude that 
even in the worst conditions this uncertainty is negligible. 

The CoG procedure is similar to the one used to compute the 
concentration index, with some minor differences. The center of 
each galaxy here is the photometric center and not the symmetry 
center. To achieve a better comparison between the two differ- 
ent bands, the Ha center was forced to coincide with the r-band 
center. Furthermore, the elliptical apertures used for Ha photom- 
etry were defined by r-band isophotes. The effective radii were 
defined as the radii containing 50% of the maximum curve of 
growth value. The effective radii were considered reliable only 
if they exceeded 1 .5 times the seeing in both bands (Arribas et al. 
2012). For AGN the central region of the Ha image was masked 
before computing of the Ha effective radius. 

4.2. EW/r parameter 

Visual inspection of the Ha images reveals that the bulge of the 
early type spirals often shows diffuse line emission. This can be 
due to ionization from AGNs, or from red and old PAGB stars 
(Stasihska et al. 2008) or a technical artifact due to improper 
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Fig. 4. Comparison between CAS parameters computed in this work and those given by Conselice (2003). The dashed lines mark 
the 1 : 1 relation. The red point is the VCC 1987 (NGC 4654). 

Table 1. Observational specifications and CAS parameters computed in this work for a set of 18 galaxies in common with Conselice 
(2003). 



NGC 


Telescope" 


Plate Scale 


Seeing 


C r 


A r 


S r 


C r 


A,. 


Sr 






(arcsec pixeL 1 ) 


(arcsec) 




This Work 


Conselice (2003) 


3596 


SPM 


0.31 


1.39 


2.79 


0.17 


0.34 


2.93 


0.14 


0.21 


4189 


SPM 


0.31 


2.20 


2.35 


0.28 


0.32 


2.46 


0.26 


0.37 


4192 


OHP 


0.69 


2.62 


3.67 


0.22 


0.64 


3.62 


0.19 


0.36 


4216 


OHP 


0.69 


3.51 


5.16 


0.21 


0.41 


5.08 


0.17 


0.38 


4254 


OHP 


0.69 


1.93 


3.24 


0.31 


0.48 


3.32 


0.29 


0.49 


4303 


CA 


0.41 


1.59 


2.86 


0.25 


0.59 


2.79 


0.24 


0.51 


4321 


OHP 


0.69 


2.10 


2.94 


0.16 


0.38 


2.91 


0.12 


0.28 


4501 


OHP 


0.69 


2.41 


3.27 


0.16 


0.21 


3.22 


0.13 


0.20 


4535 


OHP 


0.69 


2.65 


2.48 


0.12 


0.31 


2.54 


0.12 


0.40 


4548 


CA 


0.41 


1.64 


3.92 


0.07 


0.13 


3.71 


0.00 


0.03 


4569 


OHP 


0.69 


3.86 


3.78 


0.15 


0.22 


3.36 


0.08 


0.16 


4579 


CA 


0.41 


1.64 


4.18 


0.06 


0.07 


4.04 


0.03 


0.06 


4654 


CA 


0.41 


1.68 


2.73 


0.30 


0.41 


2.83 


0.11 


0.17 


4689 


CA 


0.41 


1.72 


2.74 


0.01 


0.20 


3.05 


0.03 


0.08 


5248 


SPM 


0.31 


1.79 


3.48 


0.22 


0.42 


3.46 


0.16 


0.26 


5364 


SPM 


0.31 


1.39 


3.07 


0.15 


0.38 


3.15 


0.12 


0.27 


5669 


SPM 


0.31 


1.44 


3.06 


0.19 


0.38 


3.19 


0.18 


0.31 


5701 


SPM 


0.31 


1.73 


3.72 


0.02 


-0.01 


4.11 


0.04 


0.06 



Notes. (a) SPM = San Pedro Martir 2. lm; CA = Calar Alto 3.5m; OHP = Observatoire de Haute Provence 1 .2m. 



normalization of the OFF-band images. Indeed, as pointed out 
by Spector et al. (2012), the normalization coefficient of the r- 
on the narrow-band image strongly depends on the color of the 
galaxy. We calibrated our normalization for blue galaxies (< g - 
r > — 0.4) because our sample is dominated by blue dwarfs. The 
issue arises for the bulges of massive spirals that can be up to 
1 mag redder. In these cases the underlying continuum may be 
under-subtracted (up to a factor of ~ 10%), leaving a residual 
in the Ha images that is completely unrelated to ongoing star 
formation. 

To overcome this technical problem and the additional un- 
certainties due to central AGN-like emission, we defined a new 
parameter, called EW/r, as the ratio between the effective radii 
computed in the equivalent width of the Ha (EW) map and in the 
/--band image. By dividing the Ha by the r-band continuum im- 
ages 0, we obtain a map of the EW that is defined as the strength 



of a line normalized to the underlying continuum. The signal 
of emitting HII regions with a shallow continuum is enhanced 
in this image, and on the other hand the signal in the regions 
dominated by the old stellar populations, like in the bulges, is 
quenched. 

Another important correction needs to be performed before 
dividing the Ha and r images. Since the background of both im- 
ages is zero on average with a Gaussian distribution, dividing 
the two images can produce spurious high values in pixels with 
a low r-band flux. Using the IRAF task imreplace we replaced all 
pixel values from -co to +3<T SJ t y with zero. Then, by setting the 
IRAF task imarith to place a when trying to divide by zero, we 
obtained an EW image with zero background below the sensitiv- 



3 To achieve a real EW, the pixel values should be multiplied by the 
width of the narrow band filter used for the line observations, but be- 



cause we do not require absolute calibration to compute effective radii 
(they are normalized to the maximum value of the curve of growth) we 
multiplied the pixel values by the arbitrary value of 100. This was to 
avoid systematic errors due to the poor handling of small pixel values 
by Funtools. 
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Fig. 5. Gunn r band (left), Ha (center), and EW (right) images of three galaxies spanning the full range of the EW/r parameter: from 
EW/r = 0.55 (AGC 7874: top), EW/r = 1.06 (CGCG 160-098: middle), to EW/r = 3.87 (AGC 7874: bottom). The red contours 
are drawn from the Ha map, and the black line marks the 1 arcmin scale. 



ity limifj It should be noted that setting such a sensitivity limit 
could lead to an underestimate of EW/r. Nonetheless, we veri- 
fied that lowering the limit to +2<T s ky increases EW/r by ~ 10%, 
but also increases the probability that noisy regions in the back- 
ground around the galaxy are included in the galaxy flux. Since 
the EW image translates into a map of the specific star forma- 
tion rat^l this EW/r parameter traces the position of newly born 
stars weighted on the distribution of the old stars, relative to the 
position of the old stars. It is computed as follows: 



EW/r = 



r e (EW Ha ) 
r e (Gunn r) 



(6) 



If this parameter is greater than 1, new stars are born outside 
the already constituted assembly of old stars. On the other hand, 



6 The few isolated pixels greater than the +3a- sky limit have been re- 
placed with a using a Python script. 

7 The specific star formation rate (SSFR) is the star formation rate 
normalized to the stellar mass of a galaxy. 



galaxies whose EW/ r is smaller than 1 are building up their cen- 
tral concentration. Figure |5]shows the r-band, Ha, and EW im- 
ages of three galaxies spanning the full range of the EW/ r pa- 
rameter. 



5. Results 

5.1. The CAS parameter space in the r band 

When investigating the three-dimensional space of CAS param- 
eters computed in the Gunn r band we must remove edge-on 
galaxies from our sample because dumpiness (and to a lesser 
extent also asymmetry) is affected by inclination. This is because 
the presence of a projected central dust line completely changes 
the apparent morphology of edge-on galaxies as defined in the 
CAS relative to nearly face-on objects. To avoid this projection 
effect it is recommended to work with galaxies whose axial ra- 
tio is smaller than the disk oblateness. Conselice (2003) studied 
the dependence of S as a function of the axial ratio and, fol- 
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Fig. 6. Galaxy distribution (with inclinations lower than 66.5° on the sky) in the space of CAS parameters measured in the r band. 
The solid black line is the best fit of Conselice's (2003) data and the red-orange line is the best fit of our data restricted to the 
spirals (Sa-Sd) to be consistent with the Frei et al. (1996) selection used by Conselice. Points are color-coded according to the 
morphological type: red for Sa-Sb, orange for Sc-Sd, green for Sdm-Im, and blue for the BCDs. The black contoured round points 
are the average values for each group of morphological type, and the black contoured triangles are the average values given by 
Conselice (2003) for the same group of types. Error bars mark the l<x standard deviation. 



lowing his results, we set a rather conservative limit at e < 0.4, 
where e = logio(a/b) as defined in de Vaucouleurs et al. (1991). 
This limit corresponds to galaxies with an inclination lower than 
66.5° on the sky. 

Giant spirals are clearly separated in the C r vs A, panel (see 
Fig-Eb since the concentration index decreases along the Hubble 
sequence toward later types (and with decreasing stellar mass). 
The BCDs are as concentrated and asymmetric as the Sc-Sd, al- 
though they are dwarfs. The A, vs S r panel is the only one show- 



ing a linear correlation between the parameters, as already noted 
by Conselice (2003). The best bisector fit (Isobe 1990) of our 
data is 

A r = (0.48 + 0.02) xS r + (0.05 ± 0.01). (7) 

Both the slope and the intercept of the fit are higher than what 
was found by Conselice (2003), probably because of the better 
resolution and seeing of our images compared to those of Frei 
et al. (1996). This is highlighted also by the fact that only a few 
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Fig. 7. CAS parameters in the r band (for galaxies with an inclination lower than 66.5° on the sky) as a function of Hubble type. 
The red points are the average values for each type and the error bar marks the 1 cr standard deviation. 



points lie above Conselice's best fit (black solid line in Fig. [6]). 
Finally, different Hubble types are, albeit with significant scat- 
ter, separated in the S r vs C r panel. The Sa-Sb giant spirals are 
the highest concentrated and their dumpiness is similar to that 
of the dwarf irregulars. The Sc-Sd intermediate-mass spirals are 
more clumpy but less concentrated than the Sa-Sb. The irregu- 
lars are the least concentrated but still clumpy, and the BCDs are 
as concentrated as the Sc-Sd, but are definitely not clumpy (see 
Fig0. 

Table [2] gives the average values and the standard devia- 
tions of the CAS parameters in the same bins of type as re- 
ported by Conselice (2003). The agreement with previous stud- 
ies is remarkable for spirals and for the concentration of the 
irregulars. We expect the mean values of C in our work to be 
lower than those published by Conselice (2003) because we 
made use of elliptical apertures. We verified that this difference 



Table 2. Average values and Icr rms of the CAS parameters 
computed in the Gunn r band for galaxies with an inclination 
lower than 66.5° on the sky in bins of Hubble type. 



Type 




A r 


Sr 


Sa-Sb 


3.1 ±0.5 


0.10 ±0.04 


0.10 ±0.09 


Sc-Sd 


2.8 ±0.3 


0.13 ±0.06 


0.16 ±0.14 


Sdm - Im 


2.5 ± 0.3 


0.11 ±0.05 


0.09 ±0.10 


BCD 


2.7 ± 0.3 


0.13 ±0.06 


0.08 ±0.11 



in the method increases the mean values of C by ~ 0.1 - 0.3. 
Furthermore, we note that the sample used by Conselice (2003) 
consists of spiral galaxies more massive than those considered 
here (M* > 1O 9 5 M ). If we restrict our analysis to the galax- 
ies in common with Conselice (2003) (see Table [TJ and using 
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Table 3. Average values and 1 cr rms of the CAS parameters mea- 
sured on Ha images for galaxies with an inclination lower than 
66.5° on the sky in bins of Hubble type. 



Type 




A Hn 


S 


Ha 


Sa - Sb 


2.2 + 0.9 


0.28 + 0.15 


0.65 


±0.30 


Sc-Sd 


2.2 + 0.6 


0.29 ±0.17 


0.69 


±0.43 


Sdm - Im 


2.0 ± 0.9 


0.21 ±0.16 


0.45: 


±0.35 


BCD 


2.3 ±0.5 


0.36 ±0.19 


0.33 


±0.25 



circular apertures, we obtain mean values of C of 4.0 + 0.5 and 
3.0 + 0.3 for Sa-Sb and Sc-Sd. Although the statistics is insuf- 
ficient to rule out every other possible explanation, the results 
are now remarkably consistent. Hernandez-Toledo et al. (2008) 
analyzed a sample of 539 isolated galaxies from the SDSS. The 
concentration index appears to be overestimated by ~ 10% in 
their work compared to ours because they use circular apertures 
(see Sect. 13.21 1. They also found mean values for the A and S pa- 
rameters in the r band that are consistent with those in Conselice 
(2003) and in this work for Sa-Sb galaxies. 

We did not find, in our sample that dwarf irregulars are more 
asymmetric and significantly more clumpy than spirals, as found 
by Conselice (2003). Conversely, we detected a decrement of 
these values going from spirals to the irregulars, the latter popu- 
lating the region of the CAS space at S ~ 0.1. The same result 
(S ~ 0.2) was obtained by Hernandez-Toledo et al. (2008), al- 
though they did not consider late spirals (Sc-Sd) separately from 
dwarf irregulars. A possible explanation might be a different se- 
lection and definition of irregulars. Conselice (2003) used the 
dwarf irregular sample selected by van Zee (2000). This sample 
is composed of 58 galaxies selected from the UGC (Nilson 1973) 
catalog with Mb > -18 mag vega without completeness in volume 
or surface brightness. These galaxies are on average intrinsically 
brighter than the galaxies detected in HI by ALFALFA due to 

UGC incompleteness below Mb 16 mag veg - d . Inspecting the 

images shown by van Zee (2000), we note that only a few of 
her galaxies are low surface brightness objects and not all of 
them0have been analyzed by Conselice (2003). Albeit the mea- 
surement method can play a role in this issue, we expect that 
this affects the measurements by only ~ 0.1, given the scatter in 
the right panel of Fig. [4] We conclude that this discrepancy is 
not fully understood and should be addressed in detail in future 
studies using a more complete sample of irregular galaxies. 

5.2. The CAS parameter space in the Ha line emission 

The concentration index computed from Ho- images Cn a is def- 
initely not a function of Hubble type. Only a minor decrease 
can be found from Sd to Im (see Table [3). The scatter within 
each distribution is larger than the one found for C r . Hernandez- 
Toledo et al. (2008) found that the correlations of the CAS pa- 
rameters with Hubble type are scattered in any SDSS band, but 
this scatter is greatly reduced in the r— and i bandfl suggesting 
that the structure of galaxies is better revealed in the red- and 
infrared bands. 

As an example we focus on the Sa galaxies. They are bulge- 
dominated giant spirals whose color is redder than that of dwarf 

8 The only low surface brightness galaxy whose CAS parameters 
have been computed by Conselice is UGC 10669, whose dumpiness is 
nearly zero (5 = -0.07), as one might expect from the image. 

9 The z band is found to be intrinsically noisy since the poor filter 
transmission provides images with lower S/N. 



star-forming galaxies (Gavazzi et al. 2010). The red color sug- 
gests that they are dominated by old and evolved stars. They 
show ongoing star formation activity that can sometime be quite 
relevant in an absolute quantity but their star formation rate 
per unit mass (SSFR) is very low compared to the blue dwarfs 
(Gavazzi et al. 2013a, hereafter Paper II). For this reason their 
Ha concentration index spans the whole parameter range with- 
out any dependence on the morphological classification. We 
give two examples of Sa galaxies whose Cho is very different. 
AGC7538 has a starburst nucleus (confirmed both by NED and 
by the SDSS-DR7 spectrum) and low star formation on the disk 
so the concentration index in the Ha reaches Cn a = 5.1. On the 
other hand, AGC9436 has a giant bulge without star formation 
(the nuclear activity is that of a retired galaxy, Cid Fernandes 
et al. 2011; Gavazzi et al. 2011) and the star formation activity 
takes place in a thin ring in the disk giving a Cn a = 0.3 (see 
also Figure |5]bottom). The asymmetry in Ha follows a similarly 
weak trend as the one computed in the r band (but with larger 
scatter), which supports the theory that the "flocculent asymme- 
try", as it was called by Conselice et al. (2000), is caused by 
recent star formation. For the most part, the contributions to the 
asymmetry of normal galaxies come from blue star-forming re- 
gions in the arms of the spirals. Finally, the dumpiness in Ha 
follows the same trend as a function of Hubble type previously 
found in the r band, but the scatter of points within the same bin 
of type is dramatically greater (see Figure|8). The scatter is prob- 
ably due to the stochasticity that affects the star-forming regions 
especially in the dwarfs that are populated by a small number of 
HII regions (Boselli et al. 2009, Fumagalli et al. 201 1, Weisz et 
al. 2011). 

5.3. Environmental effects on effective radii 

Fig.|9]shows the ratio between effective radii (left) and the EW/r 
parameter (right) as a function of the Hi-deficiency for all galax- 
ies. The Hi-deficiency parameter (De fm) has been defined by 
Haynes & Giovanelli (1984) as the logarithmic difference be- 
tween the HI mass observed in one object and the expected value 
in isolated and unperturbed objects of similar type and size. It is 
commonly used as a proxy for the degree of the perturbation 
that galaxies experience in dense environments. In this work we 
used an updated calibration of the Hi-deficiency parameter valid 
both for giant spirals and dwarf systems (see Paper III). It is 
immediately recognizable that environmentally perturbed galax- 
ies have a more concentrated star-formation than their old stars 
(r e (Ha) < r e (r)). This also agrees with the results of Koopmann 
et al. (2006), who found, on a smaller sample of galaxies in 
the Virgo cluster, a similar correlation between HI deficiency 
and the extension of the Ha emission. More recently Cortese et 
al. (2012) have shown a strong correlation between the HI defi- 
ciency and the extension of the star-forming disk traced by the 
UV emission. They also performed a comparison between ef- 
fective and isophotal radii showing that scatter in the correlation 
with the HI content can be greatly reduced using isophotal radii 
from radial profile fitting. Unfortunately, as discussed in Sec. 
13.11 such technique can not be applied safely to Ha images. 

It has been shown by Cayatte et al. (1994), Boselli et al. 
(2006), Fumagalli & Gavazzi (2008), Welikala et al. (2008), and 
Rose et al. (2010) that the HI removal by the ram-pressure strip- 
ping (Gunn & Gott 1972) mechanism leads to a truncation of 
the HI disk. When the HI removal is so severe that the HI disk 
shrinks inside the optical radius, there is also a truncation of the 
star-forming disk. The efficiency of this process is a function 
of the environmental conditions that each galaxy is experienc- 
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Fig. 8. CAS parameters in the Ha images (for galaxies with an inclination lower than 66.5° on the sky) as a function of Hubble type. 
The red points are the average values for each type and the error bar marks the 1 cr standard deviation. 



ing (i.e., density and temperature of the intracluster medium and 
velocity of the infalling galaxy). It has been shown both by ob- 
servations and simulations that the timescale of gas stripping is 
on the order of a few 100 Myr (e.g. Yoshida et al. 2004, 2008, 
2012; Boselli et al. 2006, 2008; Kapferer et al. 2009; Tonnesen 
& Bryan 2009, 2010, 2012; Yagi et al. 2010; Fumagalli et al. 
2011). 

It is thus recommended that any investigation on the scale 
at which galaxies produce new stars at z — should be re- 
stricted to objects that are currently not undergoing environmen- 
tal disturbances: i.e., non Hi-deficient galaxies. We set a thresh- 
old between deficient and non-deficient galaxies at defni = 0.3, 
i.e., lcr above the mean deficiency of the isolated galaxies. The 
scaling relations presented hereafter are obtained only for non- 
deficient galaxies. 



5.4. Spatial growth of galaxies 

Inspecting the ratio between Ha and r effective radii as a func- 
tion of Hubble type, we find that the recently born stars (Ha) 
generally track the r-band stellar continuum, but the average val- 
ues are slightly higher than 1 for giant spirals, consistent with 1 
for late spirals and dwarf irregulars, and they become lower than 
1 for Im and BCDs (see Fig. [10] left panel). 

By comparing the Ha effective radii with those of the stel- 
lar continuum, we have two snap-shots: a current spatially re- 
solved star formation map, and an age-weighted integral map 
of the mass built in the past. These quantities, although they do 
not directly measure the growth rate of galaxies over time, can 
be used to infer the growth size at the current epoch. Those for 
which r e (Ha) > r e (r) have a star formation that is more extended 
than their assembled stellar mass, implying that at present they 
are growing in the outer regions: inside-out growth. On the other 
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Fig. 9. Ratio between effective radii (left) and the £W/r parameter (right) as a function of the Hi-deficiency. Big red dots are the 
average values along the y-axis in different bins of HI deficiency. The lcr error bar is from bootstrap resampling. 




Hubble Type 

Fig. 10. Ratio between effective radii (left) and the EW/r param- 
eter (right) as a function of Hubble type. The big red dots mark 
the average ratios and the the 1 <x error bar marks the 1 <x standard 
deviation. 



hand, galaxies for which r e (Ha) < r e {r) are now building up the 
mass in the center. 

Since the morphological classification is intrinsically subjec- 
tive, we explored possible correlations between the ratio of the 
effective radii and the concentration index for unperturbed galax- 
ies. Fig. QT| shows the effective radii of the Ha emission and 
those of the stellar continuum in the lowest and highest quartiles 
of the concentration index as measured in the r band. Points are 
color-coded in quartiles of stellar mass and the big points are the 
average ratios for each bin of stellar mass. The average ratios 
with lcr errors from bootstrap resampling are given in Table 
All averages are consistent with a 1:1 relation except those in 
the last two bins of concentration and with high masses. That is, 
the most massive and already concentrated galaxies have a star 
formation that is more extended than their stellar mass. This is 
evidence of an inside-out assembly of disk galaxies at z ~ 0. 
This result is consistent with those by Nelson et al. (2012) at 
Z~ 1 who found an average ratio < r e (Ha) / r e (r) > = 1.3 ±0.1, 
and with those of Koopmann et al. (2006) at z ~ who found 
an average ratio of 1 . 14 + 0.07, although this latter result is more 
uncertain due to small statistics and the use of the profile fit- 
ting method for the measurements. These authors also detected 



a higher ratio for the largest galaxies, whereas the less massive 
galaxies of their sample log{M*/M ) ~ 9.5 have a mean ratio of 
one, consistently with our results highlighted in Table |4] 

The newly defined EW/ r parameter is less affected by the 
uncertainties described in Sect. [4] Computing the effective ra- 
dius in the EW image makes the central emission of giant spirals 
(which we expect to be the most contaminated by [Nil] (Paper 
I) and dust extinction issues) less dominant than the outer disk 
emission, where the metallicity is lower (Zaritsky et al. 1994; 
Skillman et al.1996) and the blue color provides a better con- 
tinuum removal (Spector et al. 2012). Using EW/r increases the 
trend of the average values for unperturbed galaxies in bins of 
Hubble type than when using r e (Ha); (see Fig. [T0l right panel). 
The average values are significantly higher than 1 for Sa-Sb spi- 
rals, then they decrease toward 1 for later spirals. We show in 
the bottom panels of Fig. QT|the EW/r parameter (i.e., the ef- 
fective radius computed in the EW image as a function of the 
effective radius in the r band). Less concentrated galaxies have a 
mean value of EW/ r of one, regardless of the stellar mass. On the 
other hand, the most concentrated galaxies have a mean value of 
EW/r that strongly depends on the stellar mass. Dwarf galaxies 
lie on the 1 : 1 relation, but with increasing mass the points devi- 
ate from equality, following a relation slightly steeper than lin- 
ear. That is, concentration of the light is the main parameter that 
influences the growth of galaxies at z ~ 0, but at fixed concentra- 
tion the stellar mass plays a second-order role. Massive galaxies 
have built their central bulges at high redshift (z > 1 Nelson et 
al. 2012) and are now growing in the outer disk. Conversely, the 
dwarfs are still growing on the scale of their previously assem- 
bled mass. 

6. Conclusions and summary 

We have investigated several structural parameters of ~ 800 late- 
type galaxies in the nearby Universe in two optical bands tracing 
different stellar populations. We used the Ha emission tht traces 
the youngest stellar population in a galaxy, and the optical r band 
that traces intermediate-age and less massive stars, which con- 
stitute the bulk of the mass in galaxies at z ~ 0. The samples are 
mainly radio (21cm) selected and biased toward gas-rich sys- 
tems. This selection effect almost completely misses the early 
types which were neglected in this work and instead focuses on 
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Fig. 11. Top panels: relation between half-light radii of the Ha emission and those of the stellar continuum radiation in the lowest 
and highest quartiles of the r-band concentration index for unperturbed galaxies. Bottom panels: the same, but using the half-light 
radii of the EWn a (EW/r parameter). Points are color-coded in quartiles of stellar mass: red for log(M t /M ) > 10.3, orange for 
9.8 < log(M,/M Q ) < 10.3, green for 9.3 < log(M»/M ) < 9.8, and blue for log(M»/M G ) < 9.3. The black-contoured big points are 
the average ratios for each bin of stellar mass. 



the most highly star-forming objects in the local Universe, cov- 
ering all late-types, from giants to the low surface brightness ir- 
regulars. The HI reservoir of these galaxies is mostly unaffected 
by environmental mechanisms, and they are the most suitable 
objects for studying the effects of secular evolution on the struc- 
tural parameters. Nonetheless, we complemented the main HI 
selection with optically selected galaxies in the densest environ- 
ments to fully understand the effect of the environment on the 
typical galaxy sizes. 

A reliable statistical analysis of the CAS properties of galax- 
ies requires large, complete, and unbiased samples. When mea- 
suring the CAS parameters in the r band, it appears to be ev- 
ident that the concentration index depends on the stellar mass 
and on the Hubble type; these variables are related because most 



massive galaxies are bulge dominated, thus most concentrated. 
Whithin this sequence BCDs represent an exception since they 
are dwarfs but they are as concentrated as intermediate-mass spi- 
rals. The asymmetry and dumpiness do not exhibit strong corre- 
lations either with stellar mass or with Hubble type. Although 
noisy, we detected an increase of these two indices from gi- 
ant Sa-Sb spirals to Sc-Sd and a subsequent decrement along 
the sequence, toward later types. This agrees with Hernandez- 
Toledo et al. (2008), although they did not consider dwarf ir- 
regulars separately from late spirals. Although even noisier, this 
pattern is recognizable also in the A and S parameters measured 
in the Ha images. This evidence supports the theory suggested 
by Conselice et al. (2000) that disregarding dynamical effects 
(e.g., tidal interactions or mergers), the asymmetry of a galaxy is 
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Table 4. Average values and 1 cr variations from bootstrap resampling for the ratio r e (Ha) / r e (r) (top) or for the ratio r e (E Wu a ) I r e (f) 
(bottom) for unperturbed galaxies in quartiles of concentration index and stellar mass. No value is given for fewer than 20 objects. 





C r < 2.4 


2.4 <C r < 2.6 


2.6 < C r < 2.9 


C r > 2.9 


log(M,/M Q ) < 9.3 


0.99 


±0.04 


0.91 ±0.04 


0.99 


±0.07 




9.3 < log(M„/M ) < 9.8 


0.94 


±0.03 


0.98 ± 0.03 


0.99 


±0.03 


1.15 ±0.09 


9.8 < log(M,/M ) < 10.3 


0.98 


±0.03 


1.01 ±0.03 


1.10 


±0.03 


1.19 ±0.07 


log(M»/M ) > 10.3 






1.04 ±0.03 


1.07 


±0.03 


1.28 ±0.07 


log(M,/M Q ) < 9.3 


1.07 


±0.04 


1.02 ±0.06 


1.06 


±0.07 




9.3 < log(M„/M ) < 9.8 


1.07 


±0.03 


1.07 ±0.04 


1.08 


±0.04 


1.31 ±0.13 


9.8 < log(M,/M ) < 10.3 


1.07 


±0.03 


1.07 ±0.04 


1.24 


±0.04 


1.45 ±0.09 


log(M»/M ) > 10.3 






1.09 ±0.04 


1.18 


±0.05 


1.54 ±0.08 



dominated by the number and the distribution of its star-forming 
regions. The dependence of the CAS parameters on the age of 
the traced stellar population agrees with previous findings re- 
ported by Taylor-Mager et al. (2007), Hernandez-Toledo et al. 
(2008), and Fossati et al. (2012). Our results for the CAS pa- 
rameters in the Ha can be compared to those of Taylor-Mager et 
al. (2007) in the UV since these two bands trace similar young 
stellar populations. We found a remarkable agreement for the 
concentration index, and also for the asymmetry and dumpiness 
that are consistent within uncertainties, although it must be taken 
into account that the statistics of the UV sample used by Taylor- 
Mager et al. (2007) is limited and thus their uncertainties are 
larger. The dependence of the CAS parameters on wavelength 
is a well-known problem for high-z morphological studies. The 
optical filters generally used on the HST probe the UV rest frame 
at z > 1 ■ There have been attempts to quantify a "morphological 
k-correction" , i.e., a conversion factor from a rest-frame UV or 
optical blue band to the well-known and calibrated optical red 
bands (Taylor-Mager et al. 2007; Conselice et al. 2008). Due to 
the similar young stellar population traced by Ha with respect 
to the UV, our CAS analysis can be used as a proxy of the mor- 
phology as seen in the ultraviolet bands. Within this framework, 
the scatter in the Ha CAS parameters (see Table [3} is so severe 
that finding a proper k-correction as a function of Hubble type is 
no trivial task. Structural properties derived from rest-frame UV 
bands should be avoided as much as possbile, favoring observa- 
tions in the IR (rest-frame optical) bands. 

Analyzing the growth in size of galaxies in the Local 
Universe, we must take the effects of local overdensities into ac- 
count. It is well-known that the environment plays a fundamen- 
tal role in the galaxy evolution. Galaxies perturbed by overdense 
environments have more concentrated recent star formation than 
their old stars. This is clearly an effect of the truncation of the HI 
disk, which is the gas reservoir for the formation of new stars. 
Although a variety of mechanisms can act together to shape the 
HI content of a galaxy and its star formation in an overdense re- 
gion (see Boselli & Gavazzi 2006), combining our results with 
those in Paper III and taking into account the evidence of a gas 
depletion acting outside-in, we identified ram-pressure stripping 
(Gunn & Gott 1972) as the main process involved. Disentangling 
the environmental effects from the evolution of isolated and cur- 
rently unperturbed galaxies, we were also able to investigate the 
size at which galaxies are growing at z ~ 0. Whether the ratio 
between effective radii in the Ha and r bands or the EW/r param- 
eter were used, we found that the concentration index is the main 
parameter that describes the current growth of galaxies. Less 
concentrated objects are now growing at the same scale as their 
already assembled mass, but for high concentrations the stellar 
mass plays a second-order role. Massive galaxies have built their 
central mass at high redshift (z > 1 Nelson et al. 2012) and are 



now growing in the outer disk. This might be expected if gas is 
accreted onto a galaxy and then cools onto the galaxy disk (see 
e.g. Brooks et al. 2009). However, it is still debated if the size 
growth of massive galaxies observed from z ~ 2 down to z ~ 
is driven by mergers, as maintained by van Dokkum et al. (2010) 
or simply by the ongoing star formation. Conversely, the dwarfs, 
regardless of their concentration, are now building up their mass 
on the scale of their old stars. Recently, Perez et al. (2013), us- 
ing spatially resolved IFU data from the CALIFA (Sanchez et al. 
2012) survey, detected a similar inside-out growth for local mas- 
sive galaxies and no gradient in the age of stellar populations for 
low-mass objects, corroborating the findings of this work. 
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Appendix A: Tables of the structural parameters 

Tables Al and A2 (available only in electronic form at the CDS) 
contain the structural parameters presented in this work for the 
Local and Coma Supercluster galaxies respectively. The data are 
listed as follows: 

- Col. 1: AGC designation, from Haynes et al. (201 1); 

- Col. 2: VCC (Binggeli et al. 1985) designation. For Coma 
replaced by the CGCG (Zwicky et al. 1968) designation; 

- Cols. 3-4: RA and DEC coordinates (J2000); 

- Cols. 5-7: r-band, Ha and EWh„ effective radii (in arcsec); 

- Cols. 8-10: concentration, asymmetry and dumpiness (CAS) 
parameters from r-band images; 
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- Cols. 11-13: concentration, asymmetry and dumpiness 
(CAS) parameters from Ha images. 

The effective radii and the CAS parameters from Ha images are 
provided only when has been possible to compute them (see Sec. 
Oandgl). 



